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Abstract. We study the evolution of twin-beam propagating inside active media 
that may be used to establish a continuous variable entangled channel between 
two distant users. In particular, we analyze how entanglement is degraded 
during propagation, and determine a threshold value for the interaction time, 
above which the state become separable, and thus useless for entanglement based 
manipulations. We explicitly calculate the fidelity for coherent state teleportation 
and show that it is larger than one half for the whole range of parameters 
preserving entanglemenent. 



The crucial ingredient of quantum information is entanglement, which is the 
essential resource for quantum computing, teleportation, and cryptographic protocols. 
Entanglement is known to to be a valuable resource for improving measurements 
spectroscopy lithography 0), interferometry [S], and tomography of quantum 

devices Most of these concepts were initially developed for discrete quantum 
variables, in particular quantum bits. Recently, however, much attention has been 
devoted to the use of continuous variables (CV), especially Gaussian state of light 
by means of linear optical circuits jxj, since CV may be easier to manipulate than 
quantum bits in order to perform various information processes [S]. In the optical 
implementation of quantum information processing, the most relevant entangled 
channel is provided by the the twin-beam state (also called two-mode squeezed 
vacuum) of two modes of radiation, which, in the Fock bases, is expressed as 



\X)) = y/l-x*Y, xV \P) ® \P) ■ 



(1) 



The twin-beam is the maximally entangled state (for a given, finit e, value of energy) 
of two modes of radiation. It can be produced either by mixing two single-modc 
squeezed vacuum (with orthogonal squeezing phases) in balanced beam splitter or, 
from the vacuum, by a nondegenerate parametric optical amplifier (NOPA). The 
evolution operator of the NOPA reads as follows U\ = exp [A (a^tf — ab)] where the 
"gain" A is proportional to the interaction-time, the nonlinear susceptibility, and the 
pump intensity. We have x = tanh A, whereas the number of photons of the twin- 
beam is N = 2 sinh 2 A = 2.t 2 /(1 - x 2 ). In view of the duality squeezing/entanglement 
via balanced beam-splitter the parameter A is sometimes referred to as the squeezing 
parameter of the twin-beam . 

The Wigner function Wo(x±, y±; X2, í/2) of a twin-beam is Gaussian, and is given 
by (in the following we will omit the arguments, and denote the Wigner function by 
W) 



W = (2na 2 + 2ira 2 _) 1 exp 



(x 1 +x 2 ) 2 {V1+V2) 2 (zi-a^) 2 (í/i 
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where the variances are given by 

a\ = l/4exp{2A} a 2 _ = l/4cxp{-2A} . (2) 

In applications such teleportation or cryptography one needs to transfer entanglement 
among distant partners, and therefore to transmit entangled states along some kind 
of channel. For optical implementation this is usually accomplished by means of 
(active) optical fibers. As a matter of fact, the propagation of twin-beam in optical 
media unavoidably lead to degradation of entanglement due to decoherence induced 
by losses and noise. In this paper, we study the evolution of twin-beam in active 
optical media, such the pair of optical fibers that may be used to transmit twin-beam, 
and analyze the separability of the evolved state as a function of the fiber parameters. 
A threshold value for the interaction time, above which the entanglement is destroyed, 
will be analytically derived. 

In jÜ] and |10j the decoherence of twin-beam due to amplitude damping has been 
studied numerically in terms of the relative entropy of entanglement and the Bell 
nonlocality factor respectively. Here we analyze the propagation of twin-beam, in 
both active and passive fibers, using the Wigner function. In this way we are able to 
analytically study the evolution as well as the entanglement properties. In particular, 
we are able to check the positivity of the partial transpose (PPT), which is a necessary 
and sufficient condition for separability for two-mode Gaussian state of light. For some 
special choices of the parameters our results may be compared, and are in agreement, 
to that of Ref. ^I] where the nonlocality Bell factor is computed from the Wigner 
function, and of Ref. ^2] where the PPT condition is applied to the transformation 
induced by a genèric four-port (amplifying or absorbing) optical device. 

In the following we first study the propagation of twin-beam in active media by 
transforming the corresponding Master equation into a Fokker-Planck equation for the 
two-mode Wigner function. Then, we apply the PPT condition to check separability 
and determine a threshold value t s for the interaction time, above which the state is 
separable and thus useless for the purposes of quantum information processing. 

The propagation of a twin-beam inside active media can be modeled as the 
coupling of each part of the twin-beam with a non zero temperature reservoir. The 
dynamics can be described in terms of the two-mode Master equation 

^=Cgt = r„(l + M a )L[a]g t + T b (í + M b )L[b]g t + T a M a L[a)] Qt + T b M b L[tf] Qt (3) 

where Qt = g(t), T a = T b = T denotes the (equal) damping rate, M a — M b = M the 
number of background thermal photons, and L[0] is the Lindblad superoperator 

L[0]g t = Og t OÏ - \o^O Qt - ±Q t OO* . 

The terms proportional to L[a] and L[b] describe the losses, whereas the terms 
proportional to L[a^] and L[b^] describe a linear phase-insensitive amplification 
process. This can be due either to fiber dynamics or to thermal hopping; in both 
cases no phase information is carried. Of course, the dynamics inside the two fibers 
are independent on each other. 

The master equation can be transformed into a Fokker-Planck equation for 
the two-mode Wigner function W(xi, y\) x%, 2/2) Using the differential representation 
of the superoperators in Eq. © the corresponding Fokker-Planck equation reads as 
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follows 

<'-»- = l (í>* , ■ j + \ {jy^xi + dy^- 

where r denotes the rescaled time t = T/jt, and the drift term 7 is given by 
1 

7 ~ 2M + 1 ' 
The solution of Eq. Q can be written as 



W T = 



9 

fií . _ 

dxí /eixí, /dt/í /dt/ 2 W {x[,y[;x' 2 ,y' 2 ) \\ G T (xj\x' j )G T (y j \y' j ) 



(4) 



(5) 



(6) 



where Wo(x\,yi;X2,y2) is the Wigner function at r = 0, and the Green functions 
G T (xj\x'j) are given by 



G T (xj\x'j) 



1 



: exp 



x 'j e 
2D 2 



D" 



1 , 



") 



(7) 



Remarkably, the diffusion coefficients D 2 remains positive for all times. Eq. (0J admits 
a stationary solution, which can be easily derived from Eq. and, independently on 
the initial state, it has the Gaussian form 



■ exp 



-2 



vl 



■ exp 



-2 



n(2M + 1) V 2M + 1J n(2M + 1) ~' r \ 2M + 1 
corresponding to the (factorized) two-mode thermal density matrix given by 

Qstat — Mm <8> Mm j 

where /ím is the density matrix of a thermal state with M thermal photons 

1 / M 



Mm 



(8) 



(9) 



(10) 



1 + M \1 + M / 

Eq. (|10|l states that after a long interaction time entanglement is totally destroyed, 
independently on the amount of entanglement initially impinged into the channel. 
This is not unexpected since the losses and the amplification noise resulting from 
the propagation in the fibers induce decoherence on the (initially entangled) beams. 
The qüestions we want to answer are the following: given an initial twin-beam \X)) 
how entanglement is degraded during the propagation ? And how long does the 
entanglement survive ? In order to answer to these qüestions, we need to solve Eq. 
10 at a genèric time r. The Wigner function W T (xi,yi;x2,y2) can be obtained by 
the convolution © , which can be easily evaluated since the initial Wigner function of 
the twin-beam is Gaussian. We have 



(27TEÍ 27rI7 



exp 



{x 1 +x 2 ) 2 (yi + y 2 ) 2 {xi~x 2 ) 2 (yi -í/2) 2 



4S 2 



4E 2 



4S 2 + 



,(H) 



— ~yT 2 



with the variances given by 

) 2 S 2 _ = e-< T a 2 _ + D 2 . (12) 

A quantum state of a bipartite system is separable if its density operator can be written 
as Q — ^2kPk a k ® Tfc, where {pk} is a probability distribution and r's and u's are 
single-system density matrices. If a state is separable the correlations between the two 
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systems are of purely classical origin. A quantum state which is not separable contains 
quantum correlations i.e. it is entangled. A necessary condition for separability is the 
positivity of the density matrix g T , obtained by partial transposition of the original 
density matrix (PPT condition) I n general, PPT has been proved to be only a 
necessary condition for separability. However, for some specific sets of states PPT is 
also a sufficient condition. These include states of 2 x 2 and 2x3 Hilbert spaces |14| . 
and Gaussian states (states with a Gaussian Wigner function) of a bipartite continuous 
variable system, e.g. the states of a two-mode radiation field ^3 EU- Our analysis 
is based on this results. In fact, the Wigner function of a twin-beam produced by 
a parametric source is Gaussian and the evolution inside active fibers preserves such 
Gaussian character. Therefore, we are able to characterize the entanglement at any 
time and to give conditions on the fiber parameters to preserve entanglement after a 
given fiber length. The PPT condition on the density matrix can be rephrased as a 
condition on the covariance matrix of the Wigner function W(x\,y-y\ x 2 , y 2 ) of the two 
modes. We have that a state is separable iff 



V 



-fi > 



j 








j 



1 

-1 



where 



V pk = (AÇ P A^) = J W(Ç) 



with AÇj = Çj - (£,•) and £ = {xx,yi,x 2 ,y 2 }. 

Given a Wigner function of the form (fTTft we have (xj) — (yj) — (AxjAyk) 
0, Vj, fc and 

(Ax 2 ) = (Ay 2 ) = l/2(E 2 l· + E 2 _) (Ax,Ax 2 ) = (A yi Ay 2 ) = l/2(E 2 l· - E 2 _) , 
such that 

/ S^ + S 2 , 

V = l 



z 2 + - s 2 _ \ 

T. 2 + + T, 2 _ 

E^ + E 2 . 

£ 2 -£ 2 _ E 2 +E 2 _ 



In order to check separability, we diagonalize V + i/40 and impose that all its 
eigenvalues are greater or equal to zero. It turns out that the state described by 
W T is separable when both the variances satisfies the condition 



E 2 > - 
+ - 4 



E 2 . > 



1 



The condition E^_ > j is weaker, and thus separability is determined by the 
condition E 2 > \. which read as follows 



exp (— 7T — 2A) H — [1 — exp (—77")] > 1 . 

7 

Given the parameters M, T and A the threshold value t s above which the state becomc 
separable is given by 

!_ 7 \ / 1 --2A- 



7 \ 1 — 7è s 



(2M + 1) log 1 



2M 



(13) 
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In terms of the unrescaled time t the threshold for separability reads as 



*' = r log ( 1 + ^~J = r log ( 1+ 2M ) ■ (14) 

In Fig. we report t s (in unit 1/r) as a function of N for different vàlues of the 
thcrmal number of photons M in the channel. The threshold t s increases with N 
from zero to a saturation value given by log[l + (2M) ]/T. This means that low 
thcrmal noise, and for an excited entangled state at the input, the entanglement is 
present on a long time scale comparable with the photon lifetime in the medium. The 
corresponding threshold L s — ct s on the interaction length is of the order of many 
Kilometers assuming a small value for M (which is a reasonable approximation at room 
temperature and optical freqüències), and an attenuation factor about O.SdB / Km for 
the optical fibers. 

A question arises whether or not twin-beam after evolution can be used for 
entanglement based protocols such teleportation. If a is the input state of a teleporting 
machine the Wigner function of the output state g (the teleported state) is given by 

W[g](x,y) = JJJJdxidyidx 2 dy 2 W[X](x + x 2 ,y + y 2l x x + x 2 , -yi - y 2 ) W[a](x 2 ,y 2 ) , 

where W[a](a;, y) is the Wigner function of the input and W[X] that of the entangled 
state providing the nonlocality for teleportation. For pure state at the input a = 
\ip)(ip\, the teleportation fidelity is given by F — (ip\g\ip) and can be calculated as 
the overlap of the Wigner functions F = n J dxdy W[a](x,y) W[g](x,y). For the 
teleportation of a coherent state \ip) = \z) supported by evoluted twin-beam of the 
form Hl·lfl the fidelity is given by 

1 1 

F = 



1 + 4E 2 . 1 + e - 2X - rt + (1 - e- rt )(2M + 1) + (1 - rfj/rj ' 

In the ideal case of entanglement provided by unperturbed twin-beam we have 
E_ = er_, and thus the known formula F = (1 + e~ 2A )~ 1 . Teleportation of a coherent 
state is a truly quantum (nonlocal) protocol if the fidelity is larger than one half • 
Remarkably, the condition F > 1/2 is equivalent to the condition for entanglement 
survival Ei > |. This means that evolved twin-beam can be used for teleportation 
in the whole range of parameters preserving entanglement, i.e. that no regimes of 
unusable entanglement are present. 

The evolution and the degradation of continuous variable entanglement of twin- 
beam in optical fibers have been studied by means of two-mode Wigner function and 
PPT condition for separability. During the propagation, entanglement survives for a 
finite interval of time. In terms of the damping rate T and the thermal noise M, the 
threshold value for separability ranges from t s — to log[l + (2M) _1 ]/r, increasing 
with the initial squeezing parameter A. For low thermal noise, and for an excited 
entangled state at the input, the entanglement is present on a long time scale compared 
with the photon lifetime in the medium. Evolved twin-beam supports coherent state 
quantum teleportation in the whole range of parameters preserving entanglement. 
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Figure 1. The threshold value í 3 (in unit 1/T) as a function of the mean 
photon number of the twin-beam N for different vàlues of the thermal number 
of photons M in the channel. From top to bottom we have the curves for 
M = 0.1,0.3,0.7,1.0. 



